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Abstract 

We study a binary dynamical process that is a representation of the voter model with two 
candidates and opinion makers. The voters are represented by nodes of a network of social contacts 
with internal states labeled 0 or 1 and nodes that are connected can influence each other. The 
network is also perturbed by opinion makers, a set of external nodes whose states are frozen in 0 
or 1 and that can influence all nodes of the network. The quantity of interest is the probability of 
finding m nodes in state 1 at time t. Here we study this process on star networks, which are simple 
representations of hubs found in complex systems, and compare the results with those obtained for 
networks that are fully connected. In both cases a transition from disordered to ordered equilibrium 
states is observed as the number of external nodes becomes small. For fully connected networks the 
probability distribution becomes uniform at the critical point. For star networks, on the other hand, 
we show that the equilibrium distribution splits in two peaks, reflecting the two possible states of 
the central node. We obtain approximate analytical solutions for the equilibrium distribution that 
clarify the role of the central node in the process. We show that the network topology also affects 
the time scale of oscillations in single realizations of the dynamics, which are much faster for the 
star network. Finally, extending the analysis to two stars we compare our results with simulations 
in simple scale-free networks. 
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I. INTRODUCTION 


Network science has provided a large body of theoretical tools to investigate complex 
systems, from physics to social sciences and biology Much work has been devoted 

to the study of networks topological properties [D El a [MHl and dynamical processes 
on networks have been shown to depend sensitively on the network structure [TTH2U] . More 
recently, the response of networks to external perturbations has also been investigated mEn- 

Eg. 

Most of the networks found in nature are scale-free, characterized by a power-law degree 
distribution and by the presence of nodes whose degree greatly exceeds the average [1] . These 
special nodes, referred to as network hubs, are crucial for the structural integrity of many 
real-world systems [26], allowing for a fault tolerance behavior against random failures [8]. 
Nevertheless, if the hubs are removed from the network by an intentional attack, the network 
might fragment into a set of isolated graphs. Thus, the presence of hubs represents at the 
same time the robustness and the ‘Achilles heel’ of scale-free networks. This property has 
been extensively studied by means of percolation theory [101 EZl EH] • In addition, network 
hubs can be detected and studied using numerous different graph measures, most of which 
express aspects of node centrality [29] . 

In this paper we study the two-states Voter Model subjected to external perturbations 
in star networks and compare the results with those obtained for fully connected systems. 
The star and fully connected topologies model two extreme scenarios, corresponding to the 
presence of a single network hub and the total absence of preferentially connected nodes, 
respectively. The perturbations represent opinion makers, who have already decided who to 
vote for and whose influence extents over the entire population. They are modeled by a set 
of external nodes whose states are hxed and that connect to all nodes of the network. The 
system exhibits a phase transition from disordered to ordered states as the external pertur¬ 
bation is decreased and can be characterized by the equilibrium probability distribution of 
hnding m nodes in a given state. We show that the shape of this distribution is very similar 
for star and fully connected networks away from the phase transition, but it shows a huger- 
print of the network topology close to the critical point. For fully connected networks the 
probability distribution is uniform at the critical point, whereas for star networks it splits 
in two peaks, rehecting the two possible states of the central node. For single realizations 
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of the dynamics and weak pertnrbations the state of the network oscillates according to the 
eqnilibrinm distribntion. We show that the time scale of these oscillations is sensitive to the 
network topology, being mnch faster for the star network. We derive approximate analytical 
solntions for the star network and extend the resnlts for mnltiple stars, which can be nsed 
as a simplihed model for a scale-free network. 

In the next two sections we describe the Voter Model with opinion makers and the 
implementation of the dynamics in a general network. Exact resnlts for the fnlly connected 


network are reviewed in section IV and in section |V] we obtain the master eqnations for star 
networks and show resnlts from nnmerical simnlations. We also generalize onr calculations 
to star networks whose center contains a group of fully connected nodes, and construct an 
approximate solution for the joint effect of two network hubs, which is further compared 
with the outcome of a scale-free network. Our conclusions are presented in section |VI[ 


II. THE VOTER MODEL 

The voter model consists of a set of individuals trying to decide in which of two candidates 
to vote for [301 El] • Their opinion can be influenced by their friends, represented by a network 
of social contacts, and by opinion makers, such as journalists or politicians, whose power 
of persuasion toward one of the candidates extents over the entire population. The opinion 
makers are modeled by external ’frozen’ nodes whose states never change and that reach all 
voters equally, acting as perturbations to the network dynamics (Fig. [^. 

The intention of a voter is quantihed by its state being 0 or 1 and the number of opinion 
makers for candidates 0 and 1 are Nq and Ni respectively. At each time step a voter is 
selected at random and its state is updated: the voter can retain its opinion or adopt the 
opinion of one of its connected neighbors, which can be a friend or an opinion maker. In the 
absence of opinion makers the population eventually reaches a consensus and the network 
stabilizes with all nodes 0 or all nodes 1, which are the only absorbing conhgurations. As 
long as opinion makers are present for both candidates the network never stabilizes, but 
it does reach a statistical equilibrium where the probability that candidate 1 has a given 
number of votes becomes independent of the time. 

This dynamical process can model other interesting systems besides as an election with 
two candidates |32l |33] , such as a population of sexually reproducing (haploid) organisms 
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[Ml l35] and herding behavior in social systems [211 [36]. It is also similar to the Glanber 
dynamics of the Ising model [23, EZ] where + iVi is analogons to the temperatnre and 
Nq — Ni to an external magnetic held. 

If the nnmber of opinion makers is zero the average time to reach consensns can be 
analytically calculated in terms of the moments of the network degree distribution [38l l39] . 
However, the presence of external perturbations complicates the dynamics and solutions 
have been obtained only for simple networks and specihc distribution of frozen nodes. In 
particular, the voter model without opinion makers was studied in regular lattices where one 
individual in the population has hxed opinion (a zealot) [in|. Analytic solutions were also 
obtained for the equilibrium distribution in fully connected networks with arbitrary number 
of opinion makers in the limit where the number of voters goes to inhnity [2S1IH] • The full 
dynamical problem with hnite number of voters was hnally solved in [2S| where it was shown 
that the solution was also a good approximation for networks of different topologies, as long 
as the number of opinion makers Nq and A^i were rescaled according to the average degree 
of the network (see also 021 ). The numbers A"o and Ni were also analytically extended to 
real numbers smaller than 1, representing weak coupling between the voters and the opinion 
makers. It was shown (see also [MlllI]) that a phase transition exists between ordered states, 
where most voters have the same opinion, to a disordered state, where approximately half 
the votes go to each candidate, as Nq and A^i go from very small to very large numbers. The 
transition occurs exactly at A"o = A^i = 1 for fully connected networks of any size. Here we 
study this phase transition in the star network. 


III. NETWORK DYNAMICS 

Consider a network with N nodes specihed by the adjacency matrix A, dehned by Aij = 1 
if nodes i and j are connected and A^j = 0 otherwise. For our purpose. An = 0 (nodes do 
not connect to themselves), and for any pair of nodes it is possible to construct a path 
connecting them. Each node has an internal state which can take only the values 0 or 1. 
The nodes are also connected to Nq external nodes whose states are hxed at 0 and to Ni 
nodes whose states are hxed at 1, as illustrated in Figjl} In order to distinguish between 
the two kinds of nodes, we call the Nq + Ni external nodes fixed and the N nodes of the 
network, whose states are variable, free. Following [25] we shall treat Nq and Ni as real 
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FIG. 1. (color online) Representation of the voter model on a network. The different colors indicate 
the internal states of the voters, which can be undecided (circles), and opinion makers (squares). 
Opinion makers can affect all voters of the network but undecided voters can only influence their 
connected neighbors. 

numbers, representing weighted coupling between the opinion makers and the voters. 

The free nodes can change their internal state according to the following dynamical rule: 
at each time step a free node is selected at random and, with probability p its state remains 
the same; with probability 1—p the node copies the state of one of its connected neighbors, 
free or hxed, also chosen at random. 

Let 

X 2 ^ 2 , ■ ■ ■ 1 1 ) • ■ ■ , v} (1) 

denote a microscopic state of the network with Xj = 0 or Xj = 1 representing the state of 
node i. There is a total of 2^ possible microscopic states and we call Pt{x) the probability 
of finding the network in the state x at time t. Since a single free node can change state per 
time step, it is useful to dehne the auxiliary state x*' which is identical to x at every node 
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except at node k, whose state is the opposite of Xk, i.e., = 1 — Xk- Explicitly, 

X^ = {Xi, X2, . . . , Xk-l, 1 - Xk, Xk+l, Xn}- (2) 

With these dehnitions, the evolution equation for the probabilities can be written as 


Pt+i{x) =pPt{x) +{l-p)j^Pt{x)Y.ZiT{xi^Xi) 

1 


(3) 


+{^-p)EU^Pt{x^)nxi^x.). 

The hrst two terms take into account the probability that the network is already in state 
X and the selected node (i) does not change its state or (ii) copies the state of a neighbor 
which is identical to its own state. The last term is the probability that the network is in 
a state differing from x by a single node, which is selected and copies the state of neighbor 
opposite to its own. 

According to the dynamical rules, the transition probabilities can be written as 

1 


T{xi -)■ Xi) = 


ki + Nq + Ni 


' N 

.t=l 


Xi -Xj\+ XiNi + (1 - Xi)iVo 


and 


r(x- -)■ Xi) = 


N 

Ajj \Xi 

-t=i 


Xj \ + (1 — xi)Ni + x\Nq 


(4) 


(5) 


ki + Nq + Ni 

where ki is the degree of the node i. Using the fact that x® = 1 — x* we hnd that the two 
transition probabilities are identical and obtain 


Pt+i(x) =pPt(x) + 


(l-p) >A [P,(x) + P,(x®)] 

N 2^ ki + No + N^ 

l=l 


N 


Aij\l -Xi-Xj\+ XiNi + (1 - Xi)NQ 

.i=i 


( 6 ) 


IV. FULLY CONNECTED NETWORKS 


For networks that are fully connected the nodes are indistinguishable and the state of the 
network is fully specihed by the number m of nodes with internal state 1 [T^ E3] . Each of 
these macroscopic states corresponds to a set of A^!/[(iV — m)!m!] degenerated microscopic 
network states. Because there are only + 1 macroscopic states equations are greatly 
simplihed. The equilibrium probability ppci'^x) of Ending the network with m nodes in state 
1 is given by the Beta-Binomial distribution [25] 

T{Ni+m) T{N + No-m) 


PFcim) = A{N,NQ,N,) 


T(N — m + 1) r(m -I- 1) ’ 


(7) 
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where 


A{N,No,N,) = 


T{N + 1) T{No + iVi) 


( 8 ) 


T{N + No + N,) r(iVi) r(iVo)- 
This expression can also be written in term of Xm = m/N. In the limit N —)■ oo, Xm becomes 
a continuous variable 0 < a; < 1 and ppc converges to the Beta distribution 


The interesting feature of the solution expressed by Eq.([^ is that for iVo = iVi = 1 it gives 
PFci'n^) = l/(-^ + 1)) rneaning that all states are equally likely, as illustrated in Figj^ 

For networks of different topologies the effect of the hxed nodes is amplihed. The probabil¬ 
ity that a free node copies a hxed node is Pi = (iVo-|-A^i)/(iVo-|-A^i-|-/ci), where ki is the degree 
of the node. For fully connected networks ki = N—1 and Ppc = {No+Ni)/{No+Ni+N—1). 
For general networks an average value Pav can be calculated by replacing ki by the average 
degree. Effective numbers of hxed nodes No^f and Ni^f can be then dehned as the values of 
Nq and Ni in Ppc for which Pav = Pfc- This leads to 


Noef = fNo, = fN,, (10) 

where / = {N — l)/kav In [2S] it was shown that Eq.([^ with the above rescaling of hxed 
nodes hts very well the probability distribution for a variety of topologies. The formula was 
tested for relatively small networks of the types random, 2-D regular lattice, Barabasi-Albert 
scale-free and small world. Similar results were obtained in the context of herding behavior 
of economic agents [211112] . 


V. STAR NETWORKS 

A. Master equation 

For a star network it is convenient to set the total number of nodes to iV -|- 1. Node 1 
is the central node and it is connected to all peripheral N nodes. The peripheral nodes, 
on the other hand, are only connected to the central node. The peripheral nodes are in¬ 
distinguishable from each other and, similar to the fully connected network, there are only 
2{N + 1) macroscopic states, characterized by having m peripheral nodes in state 1 {N + 1 
possibilities) and the central node in state 1 or 0. 
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The evolution equation for the macroscopic states can be obtained from equations ([^ if 
we define ri{m,t) and rQ{m,t) as the probabilities of having m peripheral nodes in state 1 
at time t with the central node in state 1 and 0 respectively. We obtain 


ri(m, t + 1) = ri(m, t) <p + 


(l-p) 


(N + l) [ 

+ ri(m + l,t) (1 -p) 


m{Ni + 1) + (iV — m)NQ 

(TTAdTiVo) 

(m + l)A^o 


+ 


(m + Ni] 


{N + No + m)_ 


+ ri(m - l,t) (1 - p) 
+ ro(m,t) (1 -p) 


{N + 1){1 + No + N,) 
{N — m + l)(iVi + 1) 
(iV + l)(l + iVo +W) 

(m + Ni) 


{N + l){N + No + Ni)^ 


( 11 ) 


and 


ro(m, t + l) = ro(m, t) <p + 


(l-p) 


(N + l) [ 

+ ro{m + l,t) (1 -p) 


mNi + [N — m){No + 1) {N — m + Nq) 


(1 + iVi + No) 
m + l)(A^o + 1) 


{N + No + m)_ 


+ ro(m - l,t) (l-p) 
+ ri(m,t) (1 -p) 


{N + 1){1 + No + N,) 
{N — m + l)Ni 
(iV + l)(l + iVo +W) 
{N — m + No) 


{N + 1){N + No + N, 


( 12 ) 


The first two terms in these equations take into account the probability that the network 
is in the state x at time t, and to select a node that (i) does not change its state or (ii) copies 
the state of a neighbor in its own state. The last three terms represent the probability that 
the network is in a state differing from x by a single node at time t, to select this node, and 
to copy the state of a neighbor in the opposite state. 

The probability of having m nodes in state 1 in the star network is, therefore. 


ps{m,t) =ri{m-l,t)+ro{m,t). (13) 

The results provided by equations ( |IT| and (12) agree perfectly well with numerical simula¬ 
tions. A comparison with the fully connected network is shown in Figj^ The main feature 
of these results is the different way in which the transition between ordered and disordered 
states occurs: instead of the meltdown of the Gaussian distribution observed for fully con¬ 
nected networks, the Gaussian state splits in two peaks that move toward the boundaries 
m = 0 and m = N as No and Ni are decreased. 





















FIG. 2. (color online) Equilibrium probability distribution for networks with N = 100 nodes and 
different valnes of = A^o = = 10 (green dashed), A^p = 1 (thick red), Np = 0.5 (black 

dotted), Np = 0.05 (blue), (a) fully connected network; (b) star network. 

B. Approximate solutions 


The main difficulty in solving equations ( [IT| ) and ( 12 ) is that they are coupled through the 
central node. Although we have not found exact solutions, a simple enough approximation 
can be readily obtained if the central node is momentarily considered to be hxed. If the 
central node is hxed in state 1, any peripheral node sees Nq hxed nodes in state 0 and 
Ni + 1 nodes hxed in state 1. The problem reduces to that of N independent nodes. The 
asymptotic probability that a peripheral node is in state 1 is 


Therefore, the probability that m nodes are in state 1 (the central node plus m —1 peripheral 
nodes) becomes 


pi(m) = 




N—m+l 


(16) 


N 

Similarly, hxing the central node in the state 0, the asymptotic probability that a peripheral 
node is in state 1 is 

7Vi 


1 + No + Ni’ 

and the probability that m nodes are in state 1 is 

Po(m) = - z / q ) 

\m ' 


N—m 


(16) 


(17) 
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FIG. 3. (color online) Comparison between numerical simulations (thick lines) and the approximate 


equilibrium distribution Eq.(18) (thin lines). Panel (a), for Np = 10 also shows the result for a 


fully connected network with the rescaling Eq.(lO) corresponding to / = 99/2 (red dashed curve). 


The dashed blue lines in panels (b) {Np = 1) and (c) {Np = 0.1) correspond to the approximation 
described in appendix [A| The parameters are T = 2 x 10^, p = 0.5, N = 100. For the simulations 
10^ realizations were performed. Panel (d) shows the distribution for different network sizes (from 
left to right N = 100, 250, 500, 750 and 1000) for Np = y/N /4, where the peak splits in two. 


Adding these results we obtain the approximate expression 

ps{m) 


Ni ( \ ^ ^ 

Pi{m) + Po{m), 


(18) 


Nq + Ni Nq + A^i 

where we have introduced the weights A^i/(A"o + Ni) and Nq/{Nq + Ni) of the central node 
to be in the state 1 or 0, respectively. 

In this paper we will restrict our simulations to symmetric perturbations and dehne, for 


10 
























simplicity, 


(19) 


Afp = iV„ = Ni. 

Figj^ shows a comparison between simulations and the approximate formula (18). The 
two peaks are clearly related to the two states of the central node and are reasonably well 
described by the approximation. The region between the peaks is not well represented, since 
it has important contributions from flips of the central node that have been discarded. The 
dashed blue line shows the result of a better, although ad hoc, approximation described in 
the appendix that fits the entire curve with very good precision. Figj^a) also shows the 


approximation (10) obtained via rescaling of the expression for fully connected networks, 
which works well for Np » 1. 

The approximate solutions can also be used to estimate the point where the Gaussian-like 
distribution breaks in two peaks. For large N the contributions po and pi for ps become 
Gaussians centered at Nuq and Nui with variance cx^ = NNp{l + Np)/{1 + 2Np)‘^. The 
two-peak structure appears when the distance between the two centers is of the order of the 
standard deviation. This gives Np ~ y/N and numerical calculations indicate that Nc ~ 
y/N/i. Figj^d) shows the equilibrium distribution for several values of N and Np = y/N /4. 
The transition from unimodal to bimodal distribution marks the regime where the influence 
of the central node competes with the external perturbation, modifying the equilibrium 
distribution substantially with respect to the fully connected dynamics. The two peaks 
move apart slowly as the external perturbation is decreased and are clearly separated only 
when Np 1, independently of the network size N. 

Although the equilibrium distribution of states of the star network changes smoothly 
as the perturbation is decreased, the transition in behavior is rather different from what 
is observed in the fully connected network: for Np » y/N the state is disordered, with 
approximately half the nodes in state 1 and half the nodes in state 0. The standard deviation 
is cr = y/N /2 so that a/N = 1/2V^- For Np = y/N /4, when the two peak structure appears, 
the standard deviation increases by a factor of 4 to a/N = 2/y/N. As Np decreases below 
1 and the two peaks get significantly apart, the network is most likely to be found with 
either a fraction z/i = (1 -|- Np)/{1 + 2Np) or z/q = Np/{/i + 2Np) in state 1, executing fast 
collective transitions between the two states (see next subsection). This is in contrast with 
the behavior exhibited by the fully connected network, which have either most nodes 1 or 
most nodes 0 staying in each of these states for long periods of time before moving to the 
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FIG. 4. (color online) Magnetization for a single simulation for Np = 10 (green), Np = 1 (red) 
and Np = 0.05 (blue), (a) fully connected network with N = 2 x 10^ nodes; (b) star network with 
N = 2 X 10^ peripheral nodes. Time is measured in units of network size. 

other. 

C. Dynamics and Magnetization 

In analogy with the Ising model we define the average magnetization per node as 

Af = ^ - 1 (20) 

where ni is the number of nodes in state 1, so that —1<M<+1. In order to study the 
dynamics of M we run a single simulation for each network and plot M as a function of the 
time. 

Fig. 1^ shows the results for Np = 10, 1 and 0.05 (see also Fig. [^. In these plots one 
unit of time r is a Monte Carlo step, corresponding to N steps t of the dynamics, so that 
all nodes are updated, on average, at each unit of r. For the fully connected network with 
N = 20000, FigQa), M fluctuates around zero for Np = 10 (green line). The fluctuations 
increase as the critical value is approached and for Np = 1 (red line) they take the entire 
range of M. For Np = 0.05 (blue line) the system stays a substantial amount of time 
magnetized aX M = +lorM = —1, alternating from one extreme to the other. The lower 
the values of Np the longer the times the system stays in each state for a fixed value of iV, 
and similarly for increasing N for hxed Np. 
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FIG. 5. (color online) Equilibrium probability distribution ps{m) for a star network with Nc = 2 
center nodes, N = 200 peripheral nodes and Nq = Ni = Np for Np = 1 (black) and Np = 0.3 (red). 
Thick curves show the result of simulations and thin curves the approximation given by Eq.(|21[). 
The dashed blue lines correspond to the approximation described in appendix 

For the star network, FigQb), the results show two distinct features. First, the amplitude 
of the oscillations increases smoothly as Np decreases, reflecting the position of the two 
peaks of the equilibrium distribution. For Np = 1, for instance, M oscillates in the interval 
around ±0.3. Second, the oscillations are much faster, on the scale of tens of time steps 
for Np = 0.05, as compared to the thousands of time steps of the fully connected network. 
These oscillations are clearly driven by flips of the central node, which pulls the majority of 
the peripheral nodes with it. 

The large difference in the time scales displayed in Fig. can be understood from the 
network topologies. For fully connected networks the time scale measured in number of time 
steps t is well known, given by f = N{N ± 2Np — l)/2A'"p (see |25], for instance). For small 
Np we hnd r = t/N ~ N/2Np. For the star network, on the other hand, the state of the 
peripheral nodes is controlled by the central node. If the central node is in state 0, most of 
the peripheral nodes will be in state 0 as well if Np is small. The probability that the central 
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node flips from 0 to 1 can be estimated as the probability that it copies a frozen node in 
state 1: Np/{N — 1 + 2Np) ~ Np/N. The average time for this to happen is t = N/Np or 
r = t/N = 1/Np. The two time scales differ by a factor N/2, which is consistent with the 
results shown by Fig. 


D. Generalizations 


Star networks where the center is composed not by a single node, but by a group of 
totally connected nodes can also be studied within this approximation. If the center has Nc 
nodes a stationary solution can be constructed by freezing the state of the center into m 
ones and Nc — m zeros and assigning a weight to this state according to the fully connected 
distribution pFc{fn), given by Eq.([^. Equation (18) readily generalizes to 

N 


Nr. 


p{m) ^ ^pFc{k) 


k=0 


m — k 




where 


^k = 


N, + k 


( 21 ) 


( 22 ) 


N, +No + No 

and p{k)FC is given by equation ([^ with N replaced by Nc. Fig. shows an example 
with Nc = 2 where a three peak structure is clearly visible close to the phase transition 


= 1. The approximation (21) captures well the position of the peaks, but overshoots 


their height to compensate for the lost interference between the peaks. 

As a second application we consider the joint effect of two hubs in a complex network. 
If we approximate the hubs as independent star networks with a single central node, the 
probability of Ending m nodes in the state 1 is simply given by 


p{m) = j), 

i=o 


(23) 


where we have indicated explicitly the number of nodes of each star in the distribution. For 
small values of A),, the separate distributions will have two peaks, centered at, say, mi and 


Ml —mi] m 2 and M 2 — m 2 respectively. The joint distribution given by Eq.(23) will display 


four peaks at mi + m 2 , Ai — mi + m 2 , M2 — m2 + mi and Mi +M2 — mi— m2. If the hubs are 
not independent, but coupled by only a few links, we expect this peak structure to persist. 

Figure [^a) shows the stationary distribution for two independent star networks of sizes 
Ml = 101 and M 2 = 11. The splitting of the Gaussian-like peak in two occurs at A"o = A^i ~ 
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FIG. 6. (color online) (a) Equilibrium distribution for a pair of independent star networks with 
A/i = 101 and A /2 = H nodes for Np = 10 (green dashed), A^p = 1 (thick red), Np = 0.5 (black 
dotted) and Np = 0.05 (blue), (b) Distribution for a scale-free network with 100 nodes and 
Np = 0.01. The two largest hubs (shown in red in the inset) have 16 and 13 peripheral nodes and 
are connected by their centers. 

2.6 and Nq = Ni ^ 0.8, respectively. However, because the separation between the two peaks 
of the smaller star is small, its effect is felt only at much smaller values of the perturbation, 
when the two peaks of the large hub approach the borders and the distribution becomes thin. 
Fig. I^b) shows the equilibrium distribution for a more complex network with 100 nodes 
constructed with preferential attachment. The network has two main hubs (shown in red in 
the inset) with 16 and 13 peripheral nodes, respectively. The peaks in the distribution are 
signatures of the hubs. For scale-free networks with more cycles (not shown) the presence of 
the peaks is much less conspicuous and the distribution becomes again similar to the fully 
connected case with rescaled perturbations. 


VI. CONCLUSIONS 

The voter model with opinion makers is one of the simplest dynamical systems that can 
be represented on a network. It models an election between candidates where the voters are 
influenced by their social contacts and by external factors such as journalists and politicians. 
If the number of opinion makers is zero the population is certain to reach a consensus toward 
one of the candidates independently of the structure of the network. The power of the 
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opinion makers, however, depends strongly on the average degree of the network. For a 
completely connected network the transition between nearly consensus (ordered state) and 
a tie (disordered state) takes place exactly at = iVi = 1, independently of the population 
size N. For networks with average degree kav the effect of the hxed nodes is amplihed by 
a factor / = {N — l)/kav, which can be very large for natural populations. Much above or 
much below the transition from disordered to ordered states the influence of the network 
structure is negligible and only shows up in the rescaling of the opinion makers influence, 
that is large when the network is weakly connected. This is in contrast with processes 
describing the spreading of epidemics or synchronization of oscillators, where the topology 
plays a crucial role [iiHsn]. Close to the critical point, however, the network structure can 
leave signatures in the probability distribution p(m). 

For the particular case of star networks with a single central node, the Gaussian-like 
distribution displayed by ps{m) for large values of Nq and A^i splits into two peaks centered 
at A^(A^i + l)/(iVi + + 1) and NNi/{Ni + iVo + 1) reflecting the state of the central 

node being 1 or 0. The central node controls the entire system and the distribution behaves 
approximately as a single giant node with two collective states only. For iVo = iVi = 1 
the peaks are centered at 2N/3 and N/3 respectively, which is rather different from the 
distribution of fully connected networks where ppc{Tn) = 1/{N + 1) is constant. In the 
former case the election will be won by one of the candidates with approximately 67% of the 
votes, whereas in the latter, the winner can have any number of votes with equal probability. 
For small values of Nq and Ni both star and fully connected networks are likely to be found 
in ordered states, where most nodes are in state 0 or in state 1. These states, however, are 
not stable and network oscillates between the two possibilities. We found that the average 
frequency of these oscillations are much higher for star networks than for fully connected 
ones. 

When a few weakly connected hubs are present, the effects of central nodes are still 
visible, as shown by Figj^b). However, when the system is controlled by multiple hubs, 
as in a general scale-free network, the collective behavior becomes again similar to that 
predicted by the mean held approximation and the control by ‘local leaders’ becomes much 
less relevant. In these cases Eqs. ([^ and (10) provide good approximations for the equilibrium 
probability. 

As a hnal remark we note that fully connected and stars with arbitrary number of central 
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nodes seem to be the only network topologies where a simple treatment via macroscopic 
master equations similar to is possible. Even the highly symmetric ring network 

(1-D lattice with periodic boundary conditions) does not behave as if all nodes were iden¬ 
tical, since different configuration having the same number of nodes at the state 1 give rise 
to different macroscopic states. 

M.A.M.A. and D.M.S. acknowledge hnancial support from CNPq and FAPESP. C.A.M. 
was supported by CAPES. 
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Appendix A: An ad hoc approximation for the equilibrium distribution 


The approximation (18) completely discards the fact that the state of the central node 


fluctuates and fails to describe the region between the two peaks. Here we derive a better 
approximation using phenomenological ideas. We hrst dehne 

X + Ni 


v{x) = 


1 + Ao + iVi' 


(Al) 


as the equivalent of (14) and (16) for the case where the state of the central node is in the 


average state x with 0 < a; < 1. Accordingly, we dehne 

r(A + i) 


p{x, m) = 


iy{x)^-^{l - iy{x)) 


N—m+x 


(A2) 


T{N — m + X + l)r(m — x + 1) 
as the probability of hnding m nodes in state 1, including the central the peripheral nodes 
(see eqs. (|g and (|^). If c{x) is the probability distribution that the central node is in 
state X, then 

p{m) = / dx c{x, No, Ni)p{x,m). (A3) 

Jo 

For 

/V„ /V, 

(A4) 


An Al 

c{x, No, Al) = —-+ T 7 - ~ 1) 


No + Al Aq + Al 

we recover the approximation ( |I^ . 

In order to obtain better results we need to consider smoother distributions and the 
natural functional dependence for c{x) is the continuous version of ppc^ the Beta distribution 
eq.(|^: 

Here no(Ao, Ai) and ni(Ao, Ai) measure the joint effect of the external perturbations, Aq 
and Al, and of the A — 1 peripheral nodes on the central node. Because the approximation 


with the delta functions (A4) already gives a good description of the exact distribution, no 


and 77-1 should be signihcant only close to the phase transition. The choice 

n„{N„, NO = ni{No, N,) = Ar,e-(»»+«i)/2 


(A6) 


turns out to work well for all the cases tested. 

For the case of two central nodes (see Figj^ a similar procedure can be divised. We set 

X + y + Ni 


v{x,y) = 


2 + Ao + Al 


(A7) 
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with X and y representing the states of the two central nodes. The probability that m nodes 
are in state 1 becomes 


p{x,y,m) = 


so that 


r(iv + i) 


r(iV — m + X + y + l)r(m — x — y + 1) 


u{x,y)^-^-y{l - u{x,y))^-^+^+y 


p{m) = dx dyc{x,y,No,Ni)p{x,y,m). 
Jo Jo 


(A8) 


(A9) 


The probability distribution that the two central nodes are in states x and y must reproduce 


the coefficients pFc{k) in Eq.(21). Using the analogy between Eqs. (18) and (A4) it can be 


checked that the appropriate function is 


C2{x, y, No, Ni) = c{x, Nq, Ni) c{y, Nq + I - x,Ni +x). 


(AlO) 


Indeed, using the approximation (A4) for c{x,No,Ni) we see that 


C2{x,y,no,ni) ^ (jVo+jVi)(^+jVo+jVi) ^ {A^o(A^o + 1)^(2:)^(2/) 


+NoNi[5{x)5{l -y) + 5{l- x)6{y)] + A^i(A^i + 1)5(1 - x)5(l - y)} 

(All) 


whose coefficients correspond to prcik). We remark that the integrals (A3) and (A9) might 


be difficult to evaluate numerically for very small values of Nq and Ni, since the Beta 
distribution becomes very large close to x = 0 and x = 1. In this limit, however, the 
distribution is peaked close to m = 0 and m = N and the approximation provided by the 
fully connected distribution should work well. 
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